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Abstract 

For weak dispersion and weak dissipation cases, the (l+l)-dimensional KdV-Burgers equation 
is investigated in terms of approximate symmetry reduction approach. The formal coherence of 
similarity reduction solutions and similarity reduction equations of different orders enables series 
reduction solutions. For weak dissipation case, zero-order similarity solutions satisfy the Painleve 
II, Painleve I and Jacobi elliptic function equations. For weak dispersion case, zero-order similarity 
solutions are in the form of Kummer, Airy and hyperbolic tangent functions. Higher order similarity 
solutions can be obtained by solving linear ordinary differential equations. 
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I. INTRODUCTION 



inear problems arise in many fields of science and engineering. Lie group theory 

xact analytical 



Non: 

[3] greatly simplifies many nonlinear partial differential equations. 



3,9 



was thus 



solutions are nonetheless difficult to study in general. Perturbation theory 
developed and it plays an essential role in nonlinear science, especially in finding approximate 
analytical solutions to perturbed partial differential equations. 

The integration of Lie group theory and perturbation theory yields two distinct approxi- 
mate symmetry reduction methods. The first method due to Baikov, et al [3, Q] generalizes 
symmetry group generators to perturbation forms. For the second method proposed by 



Fushchich, et al 



91], dependent variables are expanded in perturbation series and approxi- 



mate symmetry of the original equation is decomposed into exact symmetry of the system 



resulted from pertur 



parison in Refs. 
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Dation. The second method is superior to the first one from the com- 
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The well known Korteweg-de Vries-Burgers (KdV-Burgers) equation 

Ut + 6uu^ + nu^xx + vu^x = 0, (1) 

with fi and u constant coefficients, is widely used in the many physical fields especially in 
fluid dynamics. The effects of nonlinearity {6uUx), dispersion {fiUxxx) and dissipation {lyuxx) 
are incorporated in this equation which simulates the propagation of waves on an elastic 

□ n 

tube filled with a viscous fluid |X2|, and the flow of liquids containing gas bubbles [13[ and 
turbulence [lijQ, etc. 

Johnson [l2|] inspected the travelling wave solutions to the weak dissipation (u <^ I) KdV- 



Burgers equation ([T]) in the phase plane by a perturbation method and developed formal 
asymptotic expansion for the solution. Tanh function method was applied to Eq. ([1]) in the 
limit of weak dispersion {jj, <^ 1) in a perturbative way In Refs. 13], perturbation 
analysis was also applied to the perturbed KdV equations 



Vt + QvVx + Vxxx + acirjxx + ac2rjxxxx + ac3{r]r]x)x = 0, a <1, (2) 

and 

Ut + uux + Uxxx = e«M + (^Puxx, e -C 1, (3) 

respectively. 
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Eq. ([T]) can also be manipulated by means of approximate symmetry reduction approach. 
Section II and section III are devoted to applying approximate symmetry reduction approach 
to Eq. ([T]) under the case of weak dissipation {u <^ 1) and weak dispersion {fx <ti 1), 
respectively. Section IV is conclusion and discussion of the results. 



II. APPROXIMATE SYMMETRY REDUCTION APPROACH TO WEAK DISSI- 
PATION KDV-BURGERS EQUATION 

According to the perturbation theory, solutions to perturbed partial differential equations 
can be expressed as a series containing a small parameter. Specifically, we suppose that the 
weak dissipation (z^ -C 1) KdV-Burgers equation ([T]) has the solution 



oo 

U 

k=0 



fe=0 

where Uk are functions of x and t, and solve the following system 

k 

Uk,t + 6 ^ Uk-iUi^x + lJ'Uk,xxx + Uk-i,xx = 0, (A; = 0, 1, ■ ■ ■ ) (5) 

i=0 

with M_i = 0, which is obtained by inserting Eq. (j4]) into Eq. ([1]) and vanishing the 
coefficients of different powers of v. 

The next crucial step is to study symmetry reduction of the above system via the Lie 
symmetry approach 



19(1 . To that end, we construct the Lie point symmetries 

ak = Xukx + Tukt - Uk, (fc = 0, 1, ■■■), (6) 

where X, T and Uk are functions with respect to x, t and Ui, {i = 0, 1, ■ ■ ■ ). The linearized 
equations for Eq. ([5]) are: 

k 

(^k,t + 6 ^ {crk~iUi,x + Uk-i(^i,x) + fi'CTk^xxx + crk-i,xx = 0, (A; = 0, 1, ■ ■ ■ ) (7) 

with (T_i = 0. Eq. ([7]) means that Eq. is invariant under the transformations Uk 
Uk + ecTk {k = 0, 1, ■ ■ ■ ) with an infinitesimal parameter e. 

There are infinite number of equations in Eqs. and ([7]) and infinite number of argu- 
ments in X, T and Uk {k = 0, 1, ■ ■ ■). To simphfy the problem, we begin the discussion 
from finite number of equations. 



Confining the range of k to {k\k = 0, 1, 2} in Eqs. ([5]), (Q and ([7]), we see tliat X, T, 
Uq, Ui and U2 are functions witli respect to x, t, uq, ui and U2- In tliis case, tlie determining 
equations can be obtained by substituting Eq. into Eq. ([7]), eliminating mq,*; ^i,t and 
M2,t in terms of Eq. and vanishing all coefficients of different partial derivatives of Uq, Ui 
and U2. Some of the determining equations read 

T = T = T = T =n 

from which we have T = T[t). Considering this condition, we choose the simplest equations 
for X 

■^uo -^Ul -^U'Z 0) 

from which we have X = X{x,t). Considering this condition, we choose the simplest 
equations for Uq, Ui and U2 

Uo^xu2 Uqmquo f^OjMofi ^0,uo«2 ^0,«iMi ^0,niM2 ^0,«2«2 ^5 

Ul^UQUQ Ul^UQUl U\^U0U2 U\^u\U\ U\u\U2 U\u2U2 0) 

^2,uoMo ^2,no«i U2^uoU2 ^2,uiJii f^2,MiJi2 U2^U2U2 0; 

of which the solutions are 

f/o = Fi(x, t)uo + F2(X, + F^{t)U2 + F4(x, t), 

Ui = F5(x, t)Mo + Fq{x, t)ui + F7(x, if:)M2 + F^ix, t), 
U2 = Fg{x,t)uo + Fio{x,t)ui + Fii{x,t)u2 + Fi2{x,t). 
Under these relations, the determining equations are simphfied to 

Xxx = F2 = F3 = F^i = Ft = Fs = Fg = Fio = F12 = Fi^x = F/i^^x = F^^t = Fq^x = -^ii.x = 0, 

Tt = 3Xx, Xt = 6F4, Fi^t = —GF^^x, Fq^i = —GF^^x, -^ii.t = —6-^4,0;, 
Tt = Xx — Fi, Tt = 2Xx — Fi + Fq, Tt = Xx + Fu — 2Fq, Tt = 2Xx — Fq + Fn, 
which provide us with 

X = Gat + cx + xo, T = Set + to, Uq = —2cuq + a, Ui = —cui, U2 = 0, 

where a, c, Xq and to are arbitrary constants. 



Similarly, limiting the range of k to {k\k = 0, 1, 2, 3} in Eqs. ([5]), and ([7]) where 
X, T, Uq, Ui, U2 and f/3 are functions with respect to x, t, uq, ui, U2 and M3, we repeat the 
calculation as before and obtain 

X = Gat + cx + Xq, T = 3ct + Iq, Uq = —2cuq + a, Ui = —cui, U2 = 0, U3 = cu^, 

where a, c, xq and to are arbitrary constants. 

With more similar calculation considered, we see that X, T and Uk {k = 0, 1, ■ ■ ■ ) are 
formally coherent, i.e., 

X = Gat + cx + xo, T = 3ct + to, Uk = {k — 2)cuk + aSk^, {k = 0, 1, ■■■) (8) 

where a, c, xq and to are arbitrary constants. The notation 5^,0 satisfying <5o,o = 1 and 
^k,o = (A; 7^ 0) is adopted in the following text. Subsequently, solving the characteristic 
equations 

da; dt d^o dt d^i dt duk dt , , 

leads to the similarity solutions to Eq. ([5]) which can be distinguished in the following two 
subcases. 



A. Symmetry reduction of the Painleve II solutions 

When c 7^ 0, for brevity of the results, we rewrite the constants a, xq and to as ca, cxq 
and cto, respectively. Solving ^ = ^ in Eq. (Q leads to the invariant 

I(x,t) =^ = {x - Sat + Xq - Sato) {3t + toy ^ . (10) 

In the same way, we get other invariants 

Io{x,t,uo) = Po = ^{St + to)k2uo-a) (11) 

and 

hix,t,Uk) = Pk = UkiSt + to)-^^''-^\ (A; = l, 2, ■■■) (12) 

from ^ = f and ^ = § {k = 1, 2, ■■■) respectively. Viewing (A; = 0, 1, ■ ■ ■ ) as 
functions of ^, we get the similarity solutions to Eq. ([5]) 

Uk = {St + to)^^'~'^n(0 + l^Skfl, ik = 0, 1, ■ ■ ■ ) (13) 
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with the similarity variable 

^ = {x -3at + xo-3ato){M + to)-^. (14) 
Accordingly, the series reduction solution to Eq. ([1]) is derived 

oo 

^ = i + E'^'(3t + to)^('-'^Pfc(0, (15) 
and the related similarity reduction equations are 

k 

I^Pkm + 6 ^'^-^.5 - ^^^-.5 + - 2)Pk + Pk-u^ = 0, (A; = 0, 1, ■ ■ ■ ) (16) 

1=0 

with P_i = 0. When k = 0, Eq. ( JT6l) is equivalent to the Painleve II type equation. The nth 
{n > 0) similarity reduction equation is actually a third order linear ordinary differential 
equation of F„ when the previous Pq, Pi, ■ ■ ■ , Pn-i are known, since Eq. f|T6|) is just 

fiPk,^^^ + 6iPoPk,c + PkPo,c) - ^Pk,c + ik- 2)Pk = fkiO, {k=l, 2, ■■■ ) (16') 
where fk is a only function of {Pq, Pi, ■ • • , Pk-i} 

k-l 

fkiO = Ffc^jPj^^ — Pk-1,^^- 

i=l 

B. Symmetry reduction of the Painleve I solutions 

When c = and to 7^ 0, we rewrite a and xq as ato and xoto, respectively. The similarity 
solutions are 

Uk = (at + ^)4,o + Pfe(0, (A: = 0, 1, ■ ■ ■ ) (17) 
6 

with the similarity variable ^ = — x + 3at^ + xot. Hence, the series reduction solution to Eq. 
(P) is 

oo 

n = at+^ + Y,u'Pk{0, (18) 

k=0 

in which Pk{^) satisfys 

k 

fiPk,^^ + 3J2 Pk-iPi - Pk-u - a6k,oC + Ak = 0, ik = 0, 1, ■ ■ ■ ) (19) 

where P_i = and are arbitrary integral constants. 
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When k = 0, Eq. (fT9|) is equivalent to the Painleve I type equation provided that a ^ 0. 
When k = 0, a = and Aq = |/i^jr)^(m^ — 1 — m'^), Eq. f|T9|) can be solved by the Jacobi 
elliptic function, 

2 

^0 = - 2m2) + 2fim'^plcn'^{pi^ + p2, m), (20) 

where pi, p2 and m are arbitrary constants. 

When A; > 0, an equivalent form of Eq. f|T9|) is 



^Pfc,« + 6PfcPo = ^7fe(0, (A; = l, 2, ■■■) (190 

where gk{0 is a function of {Po, Pi, ■ ■ ■ , as follows 

fe-i 

1=1 

From Eq. (19'), we see that Eq. f|T9l) is a second order linear ordinary differential equation 
of Pk and can be integrated out step by step when a = 0. The results read 

Pk = PoACk + k^-' j Po,i{Bk + j Po,^gkdOdC], (21) 
with arbitrary integral constants and C^. 

III. APPROXIMATE SYMMETRY REDUCTION APPROACH TO WEAK DIS- 
PERSION KDV-BURGERS EQUATION 

We search for series reduction solutions to weak dispersion (/i ^ 1) KdV-Burgers equation 
([1]). The process is similar to the section II. A system of partial differential equations 

k 

Uk,t + Q'^Uk~iUi,x + T^Uk,xx + Uk-i,xxx = 0, {k = 0, 1, ■■■) (22) 
with M_i =0, is obtained by plugging the perturbation series solution 



oo 

U = 

k=0 



X^Afc, (23) 



into Eq. ([T]) and vanishing the coefficients of different powers of fi. 
It is easily seen that the linearized equations related to Eq. (1221) are 

k 

(Jk,t + Q^icTk-iUi^x + Uk-iCTi^x) + J^(yk,xx + Crk-l,xxx = 0, (j = 0, 1, ■■■) (24) 
i=0 



with cr_i = 0. 

The Lie point symmetries ([6]) satisfy the hnearized equations fl2^ under the approximate 
equations Restricting the range of A; to {A; | A; = 0, 1, 2} in Eqs. (122]) and (121), 

we see that X, T, Uq, Ui and U2 are functions with respect to x, t, uq, ui and U2- The 
determining equations are derived by substituting Eq. ([2D into Eq. (^^, ehminating uo^t, 
Ml t and M2,t in terms of Eq. (l22l) and vanishing coefficients of different partial derivatives of 
Uq, Ui and U2. Some of the determining equations are 

T = T = T = T =0 

from which we get T = T{t). Considering this condition, the simplest equations for X in 
the determining equations are 

^110 -^Ul ■^U2 0) 

from which we get X = X{x, t). Considering this condition, we select the simplest equations 
for Uq and Ui 

Ul^XUO Ul^UQUQ Ul^UQUl U\u2 0) 

with the solution Uq = Fi{x,t)uo + F2{x,t) and Ui = F2,{t)uo + F^{x,t)ui + F^{x,t). From 
the reduced determining equations, the simplest equations for U2 read 

HOMO U2,UQUl ^2,lt01t2 ^2,UlMl f^2,MlM2 f^2,M2M2 ^' 

leading to U2 = Fq{x, t)uo + Ft{x, t)ui + Fs{x, t)u2 + Fg{x, t). 

Combined with these conditions, the determining equations are simplified to 

Xxx = F2, = F^ = Fq = Ft = Fq = Fi^x = F2^xxx = Fi,x = Fs,x = 0, 

Tt = 2Xx, Xt = 6F2, Fi^t = —6-^2,1, F2^t = —i^F2^xx, -^4,t = —6-^2,2;, Fg^t = —QF2^x, 
Tt = Xx ~ Fi, Tt = Xx — 2F4 + Fg, Tt = 3Xx — Fi + F4, Tt = 3Xx — F4 + Fg, 
It is easily seen that 

X = Gat + cx + xq, T = 2ct + to, Uq = —cuq + a, Ui = — 2cmi, U2 = —3cu2, 

where a, c, Xq and to are arbitrary constants. 
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In like manner, we obtain 

X = Gat + cx + Xq, T = 2ct + t^, Uq = —cuq + a, Ui = —2cui, U2 = —3cu2, f/3 = —Acu^,, 

where a, c, xq and to are arbitrary constants. 

Repeating similar calculation several times, we summarize the solutions to the determin- 
ing equations 

X = Gat + cx + xo, T = 2ct + to, Uk = —{k + l)cuk + adkfi, {k = 0, 1, ■■■) (25) 

where CI, C, Xq and to are arbitrary constants. The similarity solutions to Eq. fl22p from 
solving the characteristic equations are discussed in the following two subcases. 

A. Symmetry reduction of the Kummer function solutions 

When c 7^ 0, for brevity of the results, we rewrite the constants a, xq and to as ca, cxq 
and cto, respectively. Solving ^ = ^ in Eq. (Q results in the invariant 

I(x,t) =^ = {x-6at-6ato + xo){2t + toy^. (26) 

Likewise, we get other invariants 

Ioix,t,uo) = Po = i2t + to)^2(uo-a) (27) 

and 

4(a;,t,Mfc) = Pfc = Mfc(2t + to)^^'+'^ (A: = l, 2, ■■■) (28) 

from ^ = f and ^ = f (A; = 1, 2, ■ ■ ■ ) respectively. Considering {k = 0, 1, ■ ■ ■ ) as 
functions of ^, we get the similarity solutions to Eq. (122!) 

Uk = {2t + to)"^^'"^'^/^fc(0 + a4,o, (fc = 0, 1, ■ ■ ■ ) (29) 

with the similarity variable i = {x — Gat — 6atQ + xo){2t + to)~^, and the series reduction 
solution to Eq. ([1]) is 

00 

u = a + J2 f^'C^t + to)'^^'+'^n(0, (30) 

k=0 

where Pk{0 are subject to 

k 

uPk^^^ + 6 J2 Pk^^PiA - ^Pk,i -{k+ l)Pk + Pk^im = 0, (A; = 0, 1, ■ ■ ■ ) (31) 

1=0 
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with P_i = 0. 

When A; = 0, Eq. fl3Tl) has the Kummer function solution 

p ^ (3Ci - l)//[3CiC2Ki(|(l - CO, I, g) - 2K2(|(1 - Ci), I, g)] ^ Cii. 

6aC2Ki(i(l-3C0,|,g) + K2(i(l-3Ci),i,g)] ^ ' 

where Ci and C2 are arbitrary constants, and the two types of Kummer functions Ki(p, g, 
and K2(p, g, 2;) solve the differential equation 

zy"{z) + {q- z)y\z) - py{z) = 0. 

When (A; > 0), we rearrange the terms in Eq. fl3ip as 

fc-i 

uPk,ii+6{PoPk,i+PkPo,i)-^Pk,^-{M)Pk = -Q^Pk~iPr,i,-Pk-im^ (fc = 1, 2, ■■■) (33) 

i=l 

which is a second order linear ordinary differential equation of Pk when the previous Pq? -Pi, 
• • • , Pk-i are known. 

B. Symmetry reduction of Airy function and hyperbolic tangent function solutions 

When c = and to 7^ 0, we rewrite the constants a and Xq as at^ and XqIq, respectively. 
It is easily seen that the similarity solutions are 

Uk = (at + ^)4,o + PkiOAk = 0, 1, ■■■) (34) 
o 

with the similarity variable ^ = —x + 3at^ + xot, and the series reduction solution to Eq. ([T]) 
is 

00 

u = at + ^ + Y,fi'P,{0, (35) 

fc=0 

with PkiO satisfying 

k 

uPk,^ - 3 ^ Pk^.Pi - Pk-i,^^ + aSk,o^ + Ak = 0, {k = 0, 1, ■ ■ ■ ) (36) 

i=0 

where P_i = and Ak are integral constants. 

When k = and a = —b, we get the Airy function solution to Eq. fl36l) 

p ^ (3H^[giAi(l,3i(H'^(A -bO) + Bi(l,3l(H"^(A - bp)] 
° 3[CiAi{3^biy)-l{Ao-bC)) + Bi{3^bu)-l{Ao-bO)] 
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where Ci is an arbitrary constant. The Airy wave functions Ai(2;) and Bi(2;) are hnearly 
independent solutions for y{z) in the equation y"{z) — zy{z) = 0. Ai(n, z) and Bi(n, z) are 
the nth derivatives of Ai{z) and Bi(z) evaluated at z, respectively. 
The hyperbolic tangent function solution of traveling wave form 



Po = — — ptanh[ — + d)], 

with d an arbitrary constant, can be obtained from a = and p 
When /c > 0, an equivalent form of Eq. ( |36ll is 



(38) 



(36') 



where 



k-l 



9k{0 — 3 Pk-iPi + Pk-l,^^ — ^k- 



i=l 



From Eq. (36'), it is easily seen that the kth similarity reduction equation in Eq. fl36l) is 
a first order linear ordinary differential equation of Pk- Furthermore, all equations in Eq. 
(36') can be solved step by step. The results read 



Pk = exp(- / FodO 



, (A; = l, 2, 



(39) 



where are arbitrary integral constants. 



IV. CONCLUSION AND DISCUSSION 



In summary, by applying the approximate symmetry reduction approach to (1+1)- 
dimensional KdV-Burgers equation under the condition of weak dispersion and weak dis- 
sipation, we have found that the similarity reduction solutions and similarity reduction 
equations of different orders are coincident in their forms. Therefore, we summarize the 
series reduction solutions and general formulae for the similarity equations. 

For weak dissipation case, zero-order similarity solutions are equivalent to Painleve II, 
Painleve I type and Jacobi elliptic function solutions. For weak dispersion case, zero-order 
similarity solutions are in the form of Kummer function. Airy function and hyperbolic tan- 
gent function solutions. 

A;-order similarity reduction equations are linear ordinary differential equations with re- 
spect to Pk{C)- Especially, for the period solutions (expressed by Jacobi elliptic functions) 
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with solitary waves as a special case under weak dissipation, Airy function and hyperbolic 
tangent function solutions under weak dispersion, all the higher order similarity solutions 
can be solved simply by direct integration. 

The convergence of infinite series solutions remains as a problem and expects further 
study. The approximate symmetry reduction approach can be used to search for similar 
results of other perturbed nonlinear differential equations and it is worthwhile to summarize a 
general principle for the perturbed nonlinear differential equations holding analogous results. 
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